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Four of the five expected Goldstone modes, which will be 
eaten up by gauge fields, in neutral two-flavor color super- 
conductor are actually pseudo- Goldstone modes, and their 
degenerated mass is exactly the magnitude of the color chem- 
ical potential, which is introduced to guarantee the color neu- 
trality at moderate baryon density. 

PACS numbers: 11.30.Qc, 12.39.-x, 21.65.+f 

Since the attractive interaction in the antitriplet 
quark-quark channel in quantum chromodynamics 
(QCD), the cold and dense quark matter is believed to 
favor the formation of diquark condensate and in the su- 
perconducting phase [1] . In the idealized case at asymp- 
totically high baryon density, the color superconductiv- 
ity with two massless flavors and the color-flavor-locking 
(CFL) phase with three degenerated massless quarks 
have been widely discussed from first principle QCD cal- 
culations [2]. For physical applications, one is more inter- 
ested in the moderate baryon density region which may 
be realized in compact stars and, in very optimistic cases, 
even in heavy-ion collisions. To have a stable and macro- 
scopic color superconductor, one should take into account 
the electric and color charge neutrality condition [3,4] 
which lead to a new phase, the gapless color supercon- 
ductivity [5] or the breached pairing phase [6] . The most 
probable temperature for this new phase is finite but not 
zero [7]. In two flavor case, the color neutrality can be 
satisfied by introducing a color chemical potential fig in 
the four-fermion interaction theory at moderate baryon 
density [4] , or by a dynamic generation of a condensation 
of gluon field A® in the frame of perturbative QCD at ex- 
tremely high baryon density [8] where the back ground 
gluon field (Aq) plays the role of the color chemical po- 
tential 

It is generally accepted that in the two-flavor color 
superconductor, there will be five massless Goldstone 
bosons, corresponding to the spontaneously broken color 
symmetry from SUc(3) to SUc(2)- At moderate baryon 
density, if the charge neutrality condition is not taken 
into account, there are only three Goldstone modes [9]. 
Since the degenerated mass of two diquarks of the rest 
five collective modes is proportional to the net color 
charge Q$ of the system, it is expected that [9] one can 
recover the five massless Goldstone bosons by requiring 
color neutrality. While the Goldstone bosons will finally 
be eaten up by gauge fields through Higgs mechanism 
[10], it is necessary to check whether the five expected 
Goldstone modes are really massless. In this Letter we 
will show that four of the five expected Goldstone modes 



are actually pseudo-Goldstone modes, their degenerated 
mass is exactly the magnitude of the color chemical po- 
tential (Xg which is used to guarantee the color neutrality 
of the system. 

The pseudo-Goldstone bosons were generally discussed 
in theories with spontaneously broken local symmetries 
thirty years ago [11]. The higher-order correction leads 
to spinless bosons which behavior like Goldstone bosons 
but have a small mass. 

For a neutral two-flavor color superconductor, the 
quark chemical potential matrix 



/i = diag((j, u i, fi U 2, fj, u3 , (idi, (j,d2, fida) 



(1) 



in color and flavor space can be expressed in terms of 
baryon chemical potential /i^, electrical chemical poten- 
tial /x e , and color chemical potential /xs, 

Mui = Mu2 = - 2,u e /3 + ^ 8 /3 , 
M«3 = Mb/3 - 2/i e /3 - 2^ 8 /3 , 
Hdi = Vd2 = Mb/3 + Ai e /3 + /i 8 /3 , 



Md3 = Mfc/3 + Me/3 - 2/z 8 /3 , 



(2) 



where fit, controls the baryon number density, and /i e and 
fi$ have to be introduced to ensure the electric and color 
charge neutrality. If the SUc(3) color symmetry is not 
broken, the color chemical potential /ig has to vanish, 
otherwise it would break SUc(3) explicitly. However, 
when the color symmetry is broken spontaneously by a 
color-charged diquark condensate, fj,s does not need to 
be zero [12]. 

It is well-known that a nonzero quark-antiquark con- 
densate (ipip) spontaneously breaks the chiral symmetry 
SUl(2) (2) SUr(2) of the system, and the correspond- 
ing Goldstone bosons are the three pion mesons. How- 
ever, when the current quark mass m =/= 0, the pions 
get a small mass proportional to m, due to the ex- 
plicit chiral symmetry breaking. Similarly, when /is = 0, 
a nonzero diquark condensate (4>f a ij 5 € l: > e a ° 3 ifjjf3) spon- 
taneously breaks down the color symmetry SUc{3) to 
SUc(2), and causes five Goldstone bosons corresponding 
to the broken generators T4, T5, T@, T7 and Ts, where e lJ 
and e"^ 1 are totally antisymmetric tensors in flavor and 
color space, respectively, and it is assumed that only the 
first two colors participate in the condensate, while the 
third one does not. However, in presence of a nonzero 
the SUc(3) symmetry is explicitly broken down to 
SU C (2)®U C {1) with broken generators T 4 ,T 5 ,T 6 and 
Tj. Therefore, the four Goldstone bosons corresponding 
to the broken generators T^,T^^Tq and T7 will become 
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massive, and only the Goldstone boson corresponding to 
the broken generator Tg remains massless. 

Here we have not considered the color chemical poten- 
tial /Z3, since the neutrality for the color charge Q3 is 
automatically satisfied in the current case with diquark 
condensation in the third color direction. 

To investigate quantitatively the diquarks as pscudo- 
Goldstone bosons, we have to choose a suitable model in 
describing color superconductivity at moderate baryon 
densities. It is generally accepted that the Nambu 
Jona-lasinio model (NJL) [13] applied to quarks [14] of- 
fers a simple but effective scheme to study chiral sym- 
metry restoration [15], color symmetry spontaneously 
breaking [4,5,7,9,16,17] and isospin symmetry sponta- 
neously breaking [18-20]. In the mean field approxima- 
tion to quarks and random phase approximation (RPA) 
to mesons, one can obtain the hadronic mass spectra and 
the static properties of mesons remarkably well [14], es- 
pecially the Goldstone modes corresponding to the chiral 
symmetry spontaneously breaking [15] and to the isospin 
symmetry spontaneously breaking [20]. 

The flavor SU(2) NJL model is defined through the 
Lagrangian density 

C = -0 (i^df, + njo) ip + G s (V>-0) 2 + (^*75'r-0) 2 

+ G D (^ a ; 7 V^V^) (^hV^V^) , (3) 

where Gs and Go arc coupling constants in color singlet 
channel and color anti-triplet channel, respectively, ip c = 
C^) T and V> c = tp T C are charge-conjugate spinors, C — 
i7 2 7° is the charge conjugation matrix. 

Since we focus in this Letter on the color symmetry 
spontaneously breaking and the corresponding Goldstone 
modes, to simply the notation, we consider in the fol- 
lowing only the color symmetry spontaneously breaking 
phase with nonzero diquark condensate 



A 



^Gni^^e^p) , 



(4) 



and assume that the chiral symmetry is restored in this 
phase. 

In the mean field approximation, the quarks behavior 
like quasi-particles, and the diquark condensate is con- 
trolled by the gap equation [4] 



1 - 2G D I A = 



with the function 



J (2tt)3£i 



1- f(E%)- f{EL) 



(5) 



(6) 



where the quasi-particle energies are defined as E^ = 
E^ =F Sn, E^ = ^(Ipl ±p) 2 + A 2 with the two effective 
chemical potentials p and S/j, given by p — /-ih/3 — /i e /6 + 
Ms/3, and <fyx = fi e /2, and f(x) = 1/ (e x / T + l) is the 
Fermi-Dirac distribution function. 



To consider the color and electric charge neutralities 
we need to calculate the color and electric charge densi- 
ties which can also be expressed as summations of quasi- 
particle contributions [4], 



Q.8 



d 3 p \ - 

(2tt)3 2- 



Eh 
Ei 



[l-f(E%)-f(EL)) 



+ (f(Kz) + f(Eh))\=0, 

«. = /^E[^(i-/to-/(^-)) 



+ 3(f(E e + ) - f(EL)) - e(2f(El 3 ) 

,3 

, 



El. 

2tt 2 



(7) 



with quark energies E^ = |p| ± p, E^ 3 = |p| ± fi u ^ and 
Ed3 = |p| =1= Md3 ; where the last term of Q e is the contri- 
bution from the electron gas. Note that the quarks with 
color 3 are not involved in the diquark condensation, they 
behavior like free quarks in the color superconductivity 
phase. The diquark condensate A and color and electric 
chemical potentials /ig, and fj, e as functions of tempera- 
ture T and baryon chemical potential /if, are determined 
self-consistently by the above three coupled equations. 

We now investigate diquark and meson properties at 
finite temperature and chemical potentials. In the NJL 
model, the diquark and meson modes are regarded as 
quantum fluctuations above the mean field. The meson 
modes can be calculated in the frame of RPA [14,15]. 
When the mean field quark propagator is diagonal in 
color, flavor, and Nambu-Gorkov [21] space, for instance 
the case with only chiral condensation, the summation of 
bubbles in RPA selects its specific channel by choosing 
at each stage the same proper polarization function, a 
meson mode which is determined by the pole of the cor- 
responding meson propagator is related to its own polar- 
ization function ITmm(&) [14,15] only, 



l-2G s n M M(fc) = 0. 



(8) 



However, for the quark propagator with off-diagonal el- 
ements, like the cases of 77 and 77' meson spectrum [14], 
pion superfluidity [20], and color superconductivity con- 
sidered here, we must consider carefully all possible chan- 
nels in the bubble summation in RPA. 

In the two-flavor NJL model there are four meson 
modes, the scalar meson a and the three pseudoscalar 
mesons 7r + , n , 7r_ . In the current case considering color 
symmetry spontaneously breaking in the third direction 
in color space, there are six kinds of diquarks, Di,D2, D3 
and Di, D5, D$ constructed, respectively, by colors 2 and 
3, 1 and 3, 1 and 2, 2 and 3, 1 and 3, and 1 and 2. 

Since we restrict ourselves in the color symmetry spon- 
taneously breaking phase, the mixture among the differ- 
ent channels in the bubble summation in RPA is greatly 
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reduced. The dispersion relations for the meson modes 
a and ttq are determined by their own polarization func- 
tions, 



l-2G s lWfc) = , 
l-2G s IW (fc) = 



(9) 



while the ones for the other mesons 7r+ and 7r_ and all 
the diquarks satisfy the coupled equations, 



[1 - 2G s U n+w+ (k)] [1 - 2G s U n _ w _(k)] = 
[1 - 2G D Tl DlDl (k)} [1 - 2G D Tl DlDi (k)] = 

[i - 2G D n D2D2 (k)} [i - 2G D n DlD - 2 (k)] = o 

1 - 2G D U D3D3 (k) -2G D U D3D - 3 (k) 
-2G D U D - 3D3 {k) l~2G D U DiD - 3 {k) 



dct 



0. (10) 



For the meson and diquark masses computed through 
the above dispersion relations at k$ = M 2 and k = 0, 
one needs to know the explicit expressions II(fco, 0) only. 
Performing a relatively complicated but straightforward 
calculation including Matsubara frequency summation, 
they can be written as 

n CTCT (fc ) = u nono (k ) = Ji(fco) , 

H-rr+TT+iko) = n ff _ w _(-fc ) = -Mp e - k ) , 

nr>i£>i(fco) = nD 2 _D 2 (fc ) = u DlDl (-ko) = nD 2 _D S (-fc ) 

= I A + 2{k Q +^)K 1 (k ) , 

nDsDs^o) = n_D §D5 (-fco) 

= /a + (4A:g - 8A 2 )^ 2 (fc2) + 8k K 3 (k 2 ) , 
nD 3 i5 s (fco) = n Ds2J ,(fc ) - 8A 2 K 2 (k 2 ) , (11) 



where the functions K\ , A 2 and K 3 related to the di- 
quarks are defined as 



1 1 \ f(Ets) + f(E^) - 1 



+ 



1 + l\ i(E\ I - ■ ./■(£■ ) - i 



A- 



J (27T 



e-§ (/(£?;) + /(£?!)-!) , (12) 



with 



F±(fc , |p|) = fc + Ms T ^ T #a . 
Ff(k 2 )=Ei(k 2 -±(Eif 



(13) 



Since we will not discuss the meson masses in detail in 
this Letter, the explicit expressions of the functions J\ 



and J 2 related to the mesons are not listed here. How- 
ever, the relations among the meson dispersion relations 
shown in (11) can help us to understand the meson mass 
splitting in chiral symmetry restoration phase. From the 
first equation of (11) it is clear that the masses of a and 
7T° mesons become degenerate in the color superconduc- 
tivity phase, 



M„ = M,o 



(14) 



Remember that in presence of a nonzero electric chemi- 
cal potential ju e , the chiral symmetry SUl{2) (^) SUr(2) 
is explicitly broken down to Ul{1)®Ur(1) with gener- 
ator t 3 , and the masses of the three pion mesons are not 
the same in the chiral restoration phase. This can be 
easily seen from the difference among the three pion po- 
larization functions shown in (11). Only in the case of 
yu e = 0, the relation J 2 (— k ) = Ji(fco) results in degen- 



erated meson mass M a = M n 



M. 



7T+ 



M„_. 



Considering the gap equation (5) in the color super- 
conductivity phase and the relations among the polariza- 
tion functions U DlDl ,U D2D2 ,U DlDl and Hd 2 d 2 shown 
in (11), the two mass equations for Di,D 2 ,Di and D 2 
are both simplified as 



(k 2 -v 2 )K 1 (ko)K 1 (-k ) = . 



(15) 



Obviously, one solution is fc§ = yu|, and the other is de- 
termined by 



H(k$) = K 1 (k o )K 1 {-k o )=0. 



(16) 



From the comparison with the color charge density Qs, 
one can easily prove 



K 1 (k = -n 8 ) = 2 



Qs 
A 2 ' 



(17) 



Taking into account the color charge neutrality condition 
Qs = 0, one has 



H(k 2 = ri) = o. 



(18) 



Therefore, the masses of the four diquarks, Di,D 2 ,D\ 
and D 2 , are degenerate and exactly equal to the magni- 
tude of the color chemical potential, 



M Dl = M D2 = M Di = M D - 2 = |a* s | 



(19) 



in the color superconductivity phase. The lower panel of 
Fig.l shows /is as a function of baryon chemical poten- 
tial lib at zero temperature, calculated through solving 
the gap equation (5) and the color and electrical charge 
neutrality condition (7). There are three parameters in 
the NJL model (3) in chiral limit. The momentum cutoff 
A and the coupling constant Gs can be fixed by fitting 
the pion decay constant and chiral condensate in the vac- 
uum [14,15], and the coupling constant Gd in the diquark 
channel is taken as Gd = 3Gs/4 [4] in our calculation. 
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In this case the color superconductivity phase starts at 
fib/ 3 = 330MeV. Since the above analytic discussion 
for the meson and diquark masses does not depend on 
whether one considers electrical charge neutrality or not, 
we show also in the upper panel of Fig.l the color chemi- 
cal potential fig without considering electrical charge neu- 
trality, namely we take fi e = in the calculation. We see 
that in any case the magnitude of fi$ is small. Especially, 
when we take the both charge neutralities, fi s is only a 
few MeV in a wide region. That means, the four degen- 
erated pseudo-Goldstone bosons are almost massless in 
neutral color superconductor. 
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FIG. 1. The color chemical potential fig as a function of 
baryon chemical potential fi b at zero temperature in the color 
superconductivity phase. The upper and lower panels are, 
respectively, the results without (fi e = 0) and with electrical 
charge neutrality. 

The relations among the polarization functions 
U D3D3 ,U D3Ds ,Il DsD3 and nD 3 D 3 and their explicit ex- 
pressions shown in (11) reduce the mass equation for the 
modes D3 and Dg which are linear combinations of D3 
and D§ to 

k 2 [{kl ~ 4A 2 )^ 2 (fc 2 ) - AKl{kl)] = . (20) 
One of its solution is of course 

Mn 3 (M 5s ) = , (21) 

and the other massive mode ZJg is calculated nu- 

merically through 

{kl 4A 2 )X 2 (fc 2 ) - 4i^ 2 (fc 2 ) = . (22) 



Fig. 2 shows the mass of this heavy mode as a function 
of baryon chemical potential. It is around 1100 MeV 
and even more heavy in the case without electric charge 
neutrality (dashed line), in the color superconductivity 
phase. 

1250, . , , , 1 
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FIG. 2. The mass of the heavy diquark mode Z)g (D3) , 
calculated through (22), as a function of baryon chemical po- 
tential fii, at zero temperature in the color superconductivity 
phase. The dashed and solid lines are, respectively, the results 
without (fi e — 0) and with electrical charge neutrality. 

We have investigated the quantum fluctuations in the 
neutral two- flavor color superconductivity phase in mean 
field approximation to quarks together with the ran- 
dom phase approximation to mesons and diquarks in 
the frame of NJL model. We have shown analytically 
that, there is only one massless Goldstone boson, and 
the other four expected Goldstone bosons are actually 
pseudo-Goldstone bosons, and their degenerated mass 
is exactly the magnitude of the color chemical poten- 
tial fig which is introduced in the Lagrangian density 
to guarantee the color neutrality of the system and 
breaks explicitly the color symmetry from SUc(3) to 
SUc(2) £7c(l). By self-consistently determining the 
diquark condensate, color and electric chemical poten- 
tials, the pseudo-Goldstone mass is only a few MeV, the 
same order like the current quark mass which breaks 
explicitly the chiral symmetry of QCD. It is necessary 
to note that, while the one Goldstone and four pseudo- 
Goldstone modes will finally be eaten up by gauge fields, 
the nonzero color chemical potential fig is also coupled to 
the meson and the heavy diquark modes which will then 
be reflected in the measurable properties of the system. 
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